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ABSTRACT 

We report novel, high-angular resolution interferometric measurements that 
imply the near- infrared nuclear emission in NGC 4151 is unexpectedly compact. 
We have observed the nucleus of NGC 4151 at 2.2 fim using the two 10-meter Keck 
telescopes as an interferometer and find a marginally resolved source < 0.1 pc in 
diameter. Our measurements rule out models in which a majority of the K band 
nuclear emission is produced on scales larger than this size. The interpretation 
of our measurement most consistent with other observations is that the emission 
mainly originates directly in the central accretion disk. This implies that AGN 
unification models invoking hot, optically thick dust may not be applicable to 
NGC 4151. 

Subject headings: instrumentation: interferometers - galaxies: individual(NGC 
4151) - galaxies: Seyfert 

1. Introduction 

NGC 4151 is the nearest and brightest Seyfert and it has been identified as the "archetype" 
of Seyfert 1 galaxies (Neugebauer, Graham, Soifer, & Matthews 1990). It has also been called 
"one of the most enigmatic" of galaxies (Cassidy & Raine 1996) because the origin of the 
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near-infrared nuclear emission component, which dominates galactic near-infrared emission, 
is unclear. The extent to which the observed near-infrared emission comes directly from 
the accretion disk or jet rather than from clouds of gas or dust that reprocess energy from 
the accretion disk is critical to models that attempt to unify various observational classes of 
AGN (Antonucci 1993). Over the past 30 years, arguments for both thermal and nonthermal 
(synchrotron) emission mechanisms have been proposed (see Neugebauer, Graham, Soifer, 
& Matthews (1990) and references therein). Evidence of K band emission variability on 1-2 
month timescales (Penston et al. 1974; Rieke & Lebofsky 1981; Oknyanskij et al. 1999) and 
weak (0.1-0.5%) degrees of H and K band hnear polarization (Kemp et al. 1977; Ruiz et al. 
2003) can be explained by both thermal and nonthermal models. Currently, the prevailing 
view is that the near-infrared emission is thermal and arises from the inner edge of a dust 
torus that is heated primarily by UV and soft x-ray emission from the corona of the accre- 
tion disk surrounding a massive black hole (Antonucci 1993). In this picture, the source of 
the near-infrared emission would appear spatially unresolved by the largest single aperture 
infrared telescopes and would have low degrees of linear polarization. 

Observations at 11 (Neugebauer, Graham, Soifer, & Matthews 1990; Radomski et 
al. 2003) strongly imply the presence of a <35 pc diameter source where dust grains are 
heated to around 200 K and are responsible for the mid-infrared emission. Estimates for 
the inner radius of the dust, which could correspond to the inner edge of a torus, have 
been as small as about 0.03 pc (Rieke & Lebofsky 1981), with values in the 0.2-1 pc range 
presently favored (Antonucci 1993; Pier & Krohk 1992; Evans et al. 1993; Ruiz et al. 2003). 
In some AGN unification schemes, the torus is supposed to be optically and physically thick 
(Antonucci 1993). The radius of the broad line region (BLR) is about 0.01 pc (Clavel et al. 
1990), putting the smallest estimates of the dust torus' inner radius very close to the black 
hole and accretion disk. 

The parent galaxy for the NGC 4151 Seyfert nucleus is nearly face-on (i = 21° to the 
line of site) and has a weak radio jet at P.A.= 77° E. of N. (Mundell, Wrobel, Pedlar, & 
Gallimorc 2003). While the radio feature corresponding to the UV/optical nucleus has been 
uncertain, the candidate features contain emission knots that arc unresolved on angular 
scales of a few milliarcseconds (mas) and have a surface brightnesses of a few mJy/beam 
(Mundell, Wrobel, Pedlar, & GaUimore 2003; Ulvestad, Roy, Colbert, & Wilson 1998). The 
2MASS K band flux for NGC 4151 is 8.519 ± 0.018 magnitudes, and at least 90% of the 
flux is contained in the central 1 arcsec (Kemp et al. 1977). The bolometric luminosity of 
the AGN is about 10^^ erg s~^, and a 10^ M© black hole is thought to lie at the center of the 
BLR (Ulrich 2000). Following Mundell, Wrobel, Pedlar, & GaUimore (2003), we have taken 
the distance as 13.3 Mpc (for Hq — 75 km s~^), where 1 mas corresponds to 0.065 pc (75.3 
light days) in projection. 
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Interferometric combination of the twin Keck telescopes (Colavita & Wizinowich 2002) 
provides sensitivity to 2.2 /im emission on angular scales of ~ 1 mas. Our observations 
represent a ten-fold improvement in angular resolution when compared to previous near- 
infrared measurements of AGN and make it possible to test the subparsec-scale, near-infrared 
emission models of NGC 4151. Optical/IR interferometry has traditionally been limited to 
relatively bright stellar objects. These observations represent the first measurement of an 
extragalactic source with an optical/IR interferometer. 



The Keck Interferometer combines the light from the two 10 m Keck telescopes, sepa- 
rated by 85 m with the baseline oriented 38° (E of N). Both telescopes are equipped with 
adaptive optics systems (Wizinowich et al. 2002), which are required for K band interfer- 
ometer observations. The interferometer field of view on the sky is ~50 mas in the K band, 
as set by a single-mode fiber ahead of the detector; the interferometer is not sensitive to 
emission outside of this field of view. The fringe tracker uses a 4-bin synchronous fringe 
demodulation algorithm (Vasisht et al. 2002), similar to that used at the Palomar Testbcd 
Interferometer (Colavita et al. 1999). For the data presented here, the system operated at a 
200 Hz frame rate (the 4-bin acquisition time). 

NGC 4151 was observed with the Keck Interferometer on 2003 May 20 and 21. The data 
presented here arc from the white-light channel (Acciitcr=2.18 /im and AA ~ 0.3/im). Obser- 
vations consist of a series of interleaved integrations on the source and several calibrators. 
Each integration includes 120 seconds of fringe data followed by a background measurement. 
Included with the fringe data are periodic measurements of the single-telescope fluxes. The 
data are presented as the visibility amplitude squared, normalized such that an unresolved 
object has = 1.0. On May 20, only hmited data were collected on NGC 4151 and cali- 
brators due to weather; however, the source was detected and resolved on both nights, with 
internally consistent visibility. In the results discussed below we have used data only from 
May 21. 
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Observations and Data Reduction 



2.1. 



Observations 



2.2. 



Data Reduction 



The system visibility, the instrument response to a point source, is measured with respect 
to the calibrator stars HIP58918, HD 105925, and HIP60286. The calibrator angular sizes 
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were derived by fitting photometry from SIMBAD and 2MASS. All of the cahbrators have 
estimated angular diameters of less than 0.2 mas and are unresolved by the interferometer. 
The calibrator angular size uncertainty was set to 0.1 mas. Source and calibrator data were 
corrected for biases using sky calibrations as described by Colavita (1999) and averaged into 
blocks of 5 seconds each. The calibration procedure includes a per-scan correction for the 
bias due to mismatched fluxes between the two telescopes (which includes Strehl mismatch) 
using the single-telescope flux measurements. 

The data were then calibrated for the system visibility (Bodcn ct al. 1998). The cal- 
ibrated data points for the target source arc the average of the 5 seconds blocks in each 
integration, with an uncertainty given by the quadrature of the internal scatter and the un- 
certainty in the calibrator size. In addition to the measurement error, we have estimated any 
systematic error in the calibrated visibility to be less than 0.05 (Colavita et al. 2003). This 
systematic error was summed quadratically with the measurement and calibration errors. 
The resulting average visibility is 0.84 ± 0.06 at an average projected baseline of 82.7 m at 
a P.A. = 37°E of N. The calibrated visibilities and the system visibility estimates are shown 
in Fig. 1. 



3. Discussion 

3.1. Compact Emission 

If all the light in the fringe tracker field of view (approximately the diffraction limit of 
the single telescope) came from an unresolved point source, the interferometer calibrated V"^ 
would be unity. Our measured is relatively high, implying the source of this emission 
is small. A variety of compact flux distributions could produce the measured visibility; we 
have used three simple, geometric models to estimate the angular size of the emission. The 
models correspond to a "face-on" inchnation and are constrained to fit the measured V"^. 
If the emission is distributed as a single-component Gaussian, the FWHM = 0.98 mas ± 
0.18 mas. If the the emission is distributed as a ring, the ring inner diameter = 0.96 mas ± 
0.26 with a width = 0.06 mas, where the width is set to match the total K band flux for an 
optically thick medium with a temperature of 1500 k. If the emission is distributed as a point 
source with an extended (over- resolved) component, the point-source fraction of the flux = 
0.92 ± 0.04 of the flux. Whatever the details of the emission distribution or mechanism, 
the bulk of the emission comes from a compact (<1.52 mas at the 3a level ) region. Our 
measurements rule out any scenario with a majority of the K band nuclear emission coming 
from a region with a radius in projection larger than 0.05 pc. 
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3.2. Possible Emission Mechanisms 

We consider four possible mechanisms for the K band emission. Because of the combined 
constraints imphed by our measurement and other observations, we find two of the cases to 
be unhkely. 

Synchrotron Emission: VLB A imaging of NGC 4151 at a few mas scales (Ulvestad, Roy, 
Colbert, & Wilson 1998) reveals that its parsec-scale radio jet is subluminous; it is several 
orders of magnitude fainter than the parsec-scale radio jets imaged in nearby classical radio 
galaxies. Measurements at 18 cm and 6 cm of the nuclear region reveal a relatively flat 
spectrum [v"" with a ~ —0.3) source with a peak flux density of about 1 mJy mas~^. In 
contrast, the large near-infrared flux density of ~ 10^ mJy mas"^ is inconsistent with the 
radio fluxes and implies that the nuclear near-infrared source is unrelated to the radio source 
and, by inference, unrelated to a subparsec-scale synchrotron jet. 

Steir Cluster: Explaining the compact K band emission as stellar would require ~ 3 x 10^ O 
stars in a volume of ~ 0.001 pc^. This explanation is unlikely, as it would imply a density for 
young stars ~ 100, 000 times higher than in our own galactic center (Morris & Serabyn 1996; 
Genzel et al. 2000). A stellar component is also incompatible with the measured infrared 
intensity fluctuations (Rieke & Lebofsky 1981; Quillen et al. 2000). 

Thermal Dust: If the dust is distributed as a physically thick torus, the inner radius is 
much smaller than typically estimated. The smallest possible scale for thermal dust emission 
is set by the dust sublimation radius. Using the model of Barvainis (1987), we find dust 
with a temperature of 1900 K, corresponding to a sublimation radius of 0.05 pc, is consistent 
with our measurement. Thus, if the infrared emission is from centrally heated dust, it is 
consistent with our observations only if the dust is very hot and the K band emission is 
localized near the dust sublimation radius. 

Thermal Gas: A simple, face-on, geometrically thin accretion disk model, with T oc r^"^^^ 
(Shakura & Sunyaev 1973) and with an inner temperature of 3.5 x 10^ K (consistent with 
the soft x-ray emission (Ulrich 2000)), predicts a K band magnitude of 8.4, 95% of which 
is inside a radius of ~ 0.01 pc. This model predicts an unresolved point source, consistent 
with our size measurement and with approximately the measured near-infrared flux. 

Both the thermal dust and thermal gas models are consistent with our observations. 
However, centrally heated dust, which is thought to reprocess UV/optical radiation from the 
central engine, produces a specific signature known as the "infrared bump" (Barvainis 1987; 
Sanders et al. 1989). There is no evidence for a comparable feature in NGC 4151 (Rieke & 
Lebofsky 1981; Edelson & Malkan 1986; Alonso-Herrero et al. 2001). Also, a dust torus with 
a well-defined, but small, inner radius (corresponding to the sublimation radius) would give 
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similar time delays between fluctuations in the UV/optical flux and fluctuations in the H 
and K band flux. This is inconsistent with the analysis by Oknyanskij et al. (1999). Thus 
the absence of the "infrared bump" and the presence of the H / K time delay difference both 
support the interpretation that the bulk of the infrared flux measured with the interferometer 
arises from an unresolved accretion disk. 



4. Summary and Conclusions 

We have measured the angular diameter of the K band emission in the nucleus of NGC 
4151. Our observations rule out any emission mechanism that requires a large fraction of the 
nuclear K band emission to be produced at a radius of greater than 0.05 pc from the black 
hole. Taken in the context of other observations, we interpret our measurement as evidence 
that a majority of the K band emission comes from a central accretion disk. 

We dedicate these observations to the memory of Jim Kelley, the late project man- 
ager and tireless advocate of the Keck Interferometer. We thank Robert Antonucci, Roger 
Blandford, and Julian Krohk for useful discussions. The Keck Interferometer is funded by 
the National Aeronautics and Space Administration as part of its Navigator program. Part 
of this work was performed at the Jet Propulsion Laboratory, the California Institute of 
Technology, and the Michelson Science Center, under contract with NASA. Observations 
presented were obtained at the W.M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, the University of California, and 
NASA. The authors wish to recognize and acknowledge the very significant cultural role and 
reverence that the summit of Mauna Kca has within the indigenous Hawaiian community. 
We are most fortunate to have the opportunity to conduct observations from this moun- 
tain. The Observatory was made possible by the generous financial support of the W.M. 
Keck Foundation. This work has made use of software produced by the Michelson Science 
Center at the California Institute of Technology. This work has also made use of the SIM- 
BAD database, operated at CDS, Strasbourg, Prance, and the NASA/IPAC Infrared Science 
Archive, operated by the JPL under contract with NASA. 
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Fig. 1. — The raw (bottom), calibrated (middle) and system visibility (top) for NGC4151. 
In the raw data plot, the 5 sec data averages are shown with open circles and the scan 
averages are shown with filled squares. The measured visisbility is 0.84 ± 0.064; the error 
bars in the calibrated data contain a 5% systematic component. The hour angles for the 
averaged points are -0.03, 0.33 and 0.71 hours respectively. The projected baseline length 
ranges from 81.8 to 83.6 meters and the position angle from 39.3 to 33.9 degrees (E of N). 



